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ABSTRACT 



Context. There exist few variability studies of stars in the region in the Hertzsprung-Russell diagram between the A and B-star 
pulsational instability strips. With the aid of the high precision continuous measurements of the CoRoT space satellite, low amplitudes 
are more easily detected, making a study of this neglected region worthwhile. 

Aims. We collected a small sample of B stars observed by CoRoT to determine the origin of the dilferent types of variability observed. 
Methods. We combine literature photometry and spectroscopy to measure the fundamental parameters of the stars in the sample, and 
compare asteroseismic modelling of the light curves with (differentially rotating) spotted star models. 

Results. We found strong evidence for the existence of spots and differential rotation in HD 174648, and formulated hypotheses for 
their origin. We show that the distinction between pulsations and rotational modulation is difficult to make solely based on the light 
curve, especially in slowly rotating stars. 

Key words. Stars: oscillations; Stars: variables: early-type; Stars: fundamental parameters - Stars: individual: HD 174648 



1. Introduction 

, The cool boundary of the instability strip of slowly pulsating B 
stars (SPBs hereafter, Waelkens 1991) is not well constrained 
observationally. Excitation calculations predict that it should lie 

■ around 10000- 11 000 K dPziembowski & Pamiatnvkhlll993l: 

■ iMiglio et al l l2007ah . but this is highly dependent on parameters 
I such as metallicity, core overshooting and rotation. An obser- 
. vational determination of the cool border can help to tune the 

input physics of the evolutionary models. It is, however, diffi- 
cult to constrain the exact location of the instability strip purely 
based on observations, partially because the relevant fundamen- 
' tal parameters need to be known to an accuracy that can hardly 
be reached, but also because of the required observational ef- 
fort to collect large enough samples of high precision photom- 



* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with participation of ESA's RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. 
** Based on observations made with the ESO telescopes at La Silla 
Observatory under the ESO Large Programme LP182.D-0356 
*** Based on observations made with the Mercator Telescope, operated 
on the island of La Palma by the Flemish Community, at the Spanish 
Observatorio del Roque de los Muchachos of the Institute de Astrofisica 
de Canarias. 

^ Postdoctoral Fellow of the Fund for Scientific Research, Flanders 



etry and/or time-resolved high-resolution spectroscopy to detect 
the low-amplitude gravity modes. In early attempts to map the 
variability in this area of the Hertzsprung-Russell diagram em- 
pirically, i.e ., prior to the k nowledge of the pulsation excitation 
mechanism. iBaad j d 1 989allbh performed a search for line-profile 
variability in 22 dwarfs and giants of spectral types B8-B9.5. For 
none of these targets, periodic line-profile variations were de- 
tected. Given the relatively high fraction of late B stars - com- 
pared to earlier types - the many thousands of stars observed 
in the CoRoT ( Baglin et ap 120061) exoplanet fields and in the 
Kepler dBorucki et al.ll2O10ir field ar e very p romising to confront 
with the spectroscopic null results of lBaada f 1989a,b). However, 
the effort in gathering a homogeneous set of fundamental param- 
eters to select the late-B stars among those CoRoT and Kepler 
targets is an immense undertaking. 

Another reason to focus on late B stars is to study the co- 
existence of SPBs with Bp sta rs, HgMn stars and rapidly ro- 
tating stars dBriquet et al.ll2007l) to shed light on the underlying 
origins of their differences. Although a few bright SPBs have 
been observed wit h the MOST satellite (e.g.. lAerts et al.ll2006l; 
iGruber et al.ll2009]) . high-precision space-based observations of 
late B stars with well determined fundamental parameters are 
rare. In this respect, we explore the light curves of the late B 
stars measured in the CoRoT asteroseismology channel. 
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Fig. 1. Location of the sample in the T^s - log g diagram as de- 
termined with multicolour photometry. The right grey line rep- 
resents the blue edge of the 6 Set instability strip. The light and 
dark grey regions represent SPB instability strips with Z = 0.01 
(dark grey) and Z = 0.02 (light grey) and were taken from 
iMigho et alj ( l2007bh . 



The structure of this paper is as follows. In Sect.|2l we 
present the properties of the individual stars in the sample. In 
Sect.[3] we focus on HD 174648 and argue for the existence of 
spots on its surface. In the next section, we compare the light 
curve morphology of HD 174648 to that of the other stars in the 
sample. We end with a discussion on the origins of the differ- 
ences between them. 



2. Sample description 

The stars in the CoRoT late B star sample, being secondary as- 
teroseismology targets, are not well known. For four of the stars 
the spectral subclasses are not listed in the literature. The basic 
parameters of the complete sample so far is presented in Table[T] 
and Fig.[Tl and part of the chemical compositions of stars with 
known abundances are provided in Table [2] In the framewor k 
of a preparatory programme for CoRoT ( Solano et al.l l2005l) . 
the objects HD 46179, HD 181440, HD 182198 and HD 170 935 
were analysed spectr oscopically bv lNiemczura et aL I (l2009l) and 
iLefever et al] (|2010() . To obtain a consistent determination of 
the fundamental parameters {Tss, log g, and E{B - V)) for the 
whole sample, we fitted Kurucz atmosphere models jKuruczl 
Il993h to the literature photometry for these stars, assuming the 
average extinction cur ve for the Milky Way with /?,, - 3.1 
dChiar & TielensI 2006b and solar m etallic i ty. We used the cat- 
alogues of | H0g et al.1 (2000). Cutri et al.l (l2003h. iMonet et alJ 
(2003'), Hauck & Mermilliod (1998). Mermilliod et al.' ( Il997h . 
Thompson et al.. 01978,), .Ofek (,2008,) , .Droege et al.. (.20061) , and 
van derBlieketal.l(ll996l)T and calibrated th e photometry using 

\T „f'rj„ui:„ o. r^:^^:^ ilA^nnAl\ 1 



the Vega model of'Bohlin & Gilhlandl (l2004 and the zero points 
of|Maiz-Apellaniz (2007). The best model and 90% confidence 
intervals for the parameters were determined via an interpola- 
tion of the atmosphere grid in the three parameters Teff , log g and 
E(B-V). The Introducing a scaling factor representing the angu- 
lar diameter of the star, a statistic was used to determine the 
goodness-of-fit of ~ 1 million models randomly chosen within 
the grid boundaries. Using the same statistic and all explored 
models, 90% confidence intervals on the parameters could be 



derived. We checked these results with a Nelder-Mead mini- 
mization routine and Monte-Carlo simulations. Given that the 
inclusion of UV fluxes does not always result in a better deter- 
mination of the fundamental parameters, which is partly because 
the systematic errors are larger, but also because of poor agree- 
ment between observed and predicted colours from atmosphere 
models, we did not include them. When possible, the estimated 
parameters were cross-correlated with results obtained via dif- 
ferent methods and found in the literature: the effective tempera- 
tures were also determined via Stro mgren multicolour photome - 
try using the cali bration scheme of iMoon & Dworetskvl (Il985h : 
the calibration o f|C rame 3 ([Ml was applied to Geneva photom- 
etry. Where spectra were available the reddening was estimated 
using the diffuse interstellar bands (DIBs) at 6196 A, 6613 A and 
5780 A using the calibration of lLuna et al.l (l2008l) . We compared 
these values to the reddening found via the spectral energy dis- 
tribution (SED) fit, to ensure that the effective temperature de- 
terminations are meaningful. For all stars, we found consistent 
sets of parameters, but, given the unknown systematic uncertain- 
ties, we also list the literature values besides the one we obtained 
from SED fitting. 

As a consequence of CoRoT's target selection procedure, 
the sample of stars is fairly unbiased, with the only common 
property that the stars are late B and are located near the galac- 
tic plane. At least three binaries are included, which is ex- 
pected from the hig h binarity fraction of massive stars (e.g., 
iPreibisch et al.ll2001l) . Among the targets, there are four slow ro- 
tators, one moderate rotator and three fast rotators. The different 
surface gravities suggest that the targets are spread over the main 
sequence. 

In the following, we focus on HD 174648, which is a fast 
rotator showing clear variability with a simple frequency spec- 
trum, which is interpreted as due to spots. We then compare the 
results with the light curve analysis of the other (pulsating) stars. 



3. HD 174648 

3. 1 . Fundamental parameters 

Two spectr a are available from the fibre-fed echelle spectrograph 
HERMES dRaskin et al.ll2011h on the 1.2-m Mercator telescope 
(La Palma). The spectrograph covers a range from 3800 A to 
8800 A with a resolution of ~80000. The spectra were mea- 
sured on 2010 April 18 with a S/N of -100 and -150. A first 
estimate of the rotation velocity was o btained using the Fourier 
method (ISimon-Dfaz & Herreroll2007l) applied to the profile of 
the Mg 11 4481 A line. We estimate the uncertainties via Monte 
Carlo simulations in the following way. In 10000 data sets we 
add random Gaussian noise and vary the wavelength range on 
the left and right of the line's central wavelength according to 
a uniform distribution defined on the interval [-1,+1] A with re- 
spect to the interval chosen by eye. From the first zero of the 
Fourier transform, we obtain Veq sin / - 288 + 6kms"'; the ve- 
locity connected to the sec ond zero of the Fourier tran sform is 
V2sin/ = 156 ± 3kms~'. iReiners & Schmiti (l2002bl) derived 
a relation between the ratio of the two first zeros in the 

Fourier transform and the differential rotation rate. For slowly 
and rigidly rotating massive stars, and assuming a linear limb 
darkening law, this ratio is approximately 1.76, which is on 
the low side of our estimate ^2/91 = 1-84 ± 0.05. This sug- 
gests the presence of differential rotation with polar accelera- 
tion. However, the Mgii line is actually a narrow blend, typi- 
cally resulting in an overestimation of the Veq sin ; by 10km s"' 
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Table 1. CoRoT cool B star sample description for the asteroseismology programme of the mission. 
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Table 2. Element abundances for the stars in the CoRoT cool B star sample 








Star 


He 


C 


O Mg 


Al Si 


P S 


Ca Ti V 


Cr Fe 


Ni Ba 


HD 170935 


10.79 




8.18 


7.36 




4.64 


5.66 7.57 




HD 181440 


10.90 




7.80 


7.24 




4.36 


5.09 6.85 




HD 182198 


10.93 


8.38 


7.61 7.48 


6.56 7.34 


5.28 7.06 


6.15 4.61 4.55 


5.65 7.27 


6.45 2.56 


HD 46179 


11.04 




8.85 7.58 


7.60 




4.61 


5.79 7.29 




Sun 


10.93 


8.39 


8.66 7.53 


6.37 7.51 


5.36 7.14 


6.31 4.90 4.00 


5.64 7.45 


6.23 2.17 



Notes. The abundan ce of an element X is given as loge(A'), where loge(H) = 12 (taken from iNiemczura et al.ll2009t) . Solar abundances from 
lAsplund et al] ( l2005h . 



at these rotation velocities, and raising the qilq\ value apprecia- 
bly. We checked if we could improve the shape of the rotational 
kernel by taking more U nes into account (e.g., as was done in 
iReiners & Roverir2004ah . but this was not case: there are far too 
few lines, and the continuum estimation for these few lines is not 
reliable enough. 

To determine the effective temperature and surface gravity 
of the star, we compared the obser ved spectral Unes w ith an 
ATLAS 12 grid of synthetic spectra (iPalacios et al.ll2010h . after 
broadening them according to a Veq sin / = 288 km s"' rotational 
profile and an instrumental profile with a full width at half max- 
imum of 0.08 A. We assumed a parabolic shape of the contin- 
uum, via a Nelder-Mead Simplex minimization method, where 
we minimized the residuals of the normalised observed spectral 
profile and the trial synthetic profile, via a;(f^-like statistic. 



X' 



1 ^„ 



aA- +bA + c 



The high rotation rate smears out the large majority of the lines, 
making abundance determination impossible. The most promi- 
nent lines in the observed spectrum are the Balmer lines, the 
Mgii4481A line and the Caii3933A line. The equivalent width 
of the latter is the most sensitive one to temperature, based on 
the synthetic spectra. Moreover, the low reddening minimizes 
interstellar contributions. Therefore, we perform an independent 
check to derive the temperature only using the equivalent width 
of this feature. We calculated the equivalent width of the ob- 
served line by integrating the feature between two points outside 
the line region. The equivalent width distribution is again recov- 
ered via Monte Carlo simulations. The Icr eiToron the equivalent 
width translates to eff'ective temperatures between 10000 K and 
10750 K. Taking 3cr extends this interval to 9 750 - 13 750 K, 



but at Teff > 1 1 000 K, the Ca n feature becomes less sensitive 
to the temperature. In Fig.|2] we show the fit to a collection of 
spectral lines, with synthetic spectra corresponding to 9500 K, 
10500 K and 11500K. We conclude that different methods give 
a consistent estimate of Teff ~ 10500 K, which is also consistent 
with the appearance of a weak Hel line at 4471 A and normal 
He abundance. Finally, from the width of the Balmer wings, we 
deduce a surface gravity of log g - 4.2 + 0.3. 

On the basis of the CoRoT light curve morphology (see, e.g., 

IStra ssmeier & Bopp 1992) and the star's location in the HR di- 
agram, we investigated two processes that could possibly lie at 
the origin of cyclical modulation: pulsations and spots. 



3.2. Nonlinear pulsations 

Before performing a full frequency analysis of the CoRoT Ught 
curve of HD 174648, we corrected it for discontinuities, and 
mini mized the influence o f the satellite's orbital signal (see, 
e.g., iDegroote et al.ll20Tol) . We performed a frequency analy- 
sis with nonlinear optimization after every prewhitening stage, 
using only the last 5 frequencies found at that stage, to account 
for the closely spaced peaks. Although the full frequency analy- 
sis resulted in the detection of 28 frequencies, we only list those 
with a S/N > 4 or possible combinations in Table[3] The first two 
frequencies fi and /2 are closely spaced (A/ x 1.2/T), and two 
and one of their harmonics are recovered, respectively {fj, = 2/i , 
/4 = 2/2, and fd - 3/1). The phase diagrams of these frequen- 
cies, after prewhitening all other detections, are shown in Fig. [3] 
These, and the presence of the harmonics of these two frequen- 
cies found in the Fourier transforms, reveal unusual behaviour 
compared to the sinusoidal pha se diagrams of confirmed SPBs 
(see, e.g., |De Cat & Aertsll2002h . We cannot make any firm con- 
clusions on the two low and independent frequencies and fi. 
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Fig. 2. Comparison between model template spectra (black), synthetic spectra (red) additionally deformed according to model 
parameters resulting from the spot model fit (see text), and the observed profiles (grey) (from top to bottom and left to right: 
Can 3933 A with and He, H6, Mg ii4481A and Ha). In each panel, the top black spectrum shows the comparison to a template 
with parameters Tefr = 11 500 K and \ogg - 4.0, the middle to the template T^ff - 10500/ir and \ogg = 4.0 and the bottom to 
Teff = 9 500 K and log g - 4.0. The top and bottom red dashed lines show models (with log g - 4.0) compatible with the photometry 
with Jeff = 1 1 580 K and T^s = 10 357 K, respectively. 



because their period is long compared to the total time of ob- 
servations, and they could thus be influenced by instrumental 
effects. 

Although SPB models have eigenfrequencies in the region 
of interest, the corresponding /"-values of 5 and higher result 
in strong geometric cancellation, making high amplitudes of 
~0.l% improbable. In fast rotators, geometric cancellation is 
even more eff'ective than for slow rotators because of two addi- 
tional eflFects: first, the modes are more concentrated around the 
equator (e.g., iReese et al ] l2009h . and second, there is less flux 
coming from th ese regions d ue to gravity darkening (e.g., 5 times 
less at 0.9vciit, lAerts et aD l2004). To explain the dominant fre- 
quencies with a beating phenomenon due to pulsations, it is thus 
necessary to shift low-degree model frequencies of slowly ro- 
tating SPB stars to higher frequencies; cf. the detection of three 
pulsation modes aroun d 2.5 d~' for the rapidly rotating SPB star 
HD 121 190 (.Aerts & Kolenberg|2005h . 

Rotationally shifting the low-degree frequencies from ~1 d ' 
to ~4d~' requires a high rotational angular velocity: in the case 
of an ^ = \,m - 1 gravity mode, and only retaining first order 



effects to obtain an upper limit, we get 



The maximum possible angular break-up velocities for cool B 
stars are reached at the ZAMS (when they are most compact), 
where values up to ~3.5d ' can be expected. As a star evolves 
the angular break-up velocity decreases. The required angular 
velocity to shift the frequencies can be lowered by increasing 
the mode degree, and more importantly the azimuthal order At 
the same time, the observed frequency spectrum should lose all 
its structure, since the separation between the modes becomes 
higher This is in sharp contrast to the observed frequency spec- 
trum with the two closely spaced frequencies. 

We conclude that it is highly unlikely that the observed fre- 
quency spectrum of HD 174648 is caused by stellar pulsation: 
the effective temperature is far too low to excite pressure modes 
in that region, and the frequency spectrum is too structured to 
originate from rotationally shifted nonradial gravity modes. We 
thus investigate the possibility that the observed variability does 
not arise from pulsations, but from rotational modulation. 
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Table 3. Frequency analysis of the CoRoT light curve of HD 174648 after correcting for discontinuities and removal of the satellite's 
orbital signal. 



ID 


f 




a 
ppm 




<P 

rad/27r 


S/N 


notes 


/i 


3.79144 ± 


0.00008 


1829.3 ± 


6.5 


-0.273 + 0.003 


27.42 




/2 


3.84083 ± 


0.00007 


1356.7 ± 


3.8 


-0.313 + 0.003 


31.27 




h 


7.58282 ± 


0.00008 


627.0 ± 


2.2 


0.005 ± 0.003 


31.66 


2/i - 0.00007 


/4 


7.68346 ± 


0.00032 


156.0 ± 


2.1 


-0.050 ±0.013 


23.76 


2/2-0.0018 


fi 


0.23672 ± 


0.00112 


564 ± 


2.1 


0.296 ± 0.047 


5.88 




fe 


11.37273 ± 


0.00116 


41.5 ± 


2.1 


-0.029 ± 0.050 


10.64 


3/i - 0.0016 


fi 


0.30554 ± 


0.00131 


43.1 ± 


2.1 


-0.309 ± 0.055 


5.29 




h 


3.95753 ± 


0.00153 


33.0 ± 


2.1 


-0.003 ± 0.065 


5.80 




h 


0.09974 ± 


0.00179 


26.3 ± 


2.0 


-0.021 ±0.075 


4.09 




fn 


7.63157 ± 


0.00205 


25.0 ± 


2.0 


-0.055 ± 0.086 


6.26 


/■1 + /2- 0.0007 


fn 


0.54299 ± 


0.00229 


19.5 ± 


2.0 


0.192 ±0.097 


3.65 


/,+/7 - 0.0007 


/l6 


0.27052 ± 


0.00303 


21.0 ± 


2.0 


0.110 ±0.128 


3.19 


/i3 - fis - 0.0003 


fit, 


6.02177 ± 


0.00345 


144 ± 


2.0 


-0.786 ±0.145 


3.82 


2/17 


fl9 


12.81640 ± 


0.00354 


13.5 ± 


2.0 


0.530 ±0.149 


4.13 





Notes. The frequencies, amplitudes and phases were optimized via nonlinear fitting after every prewhitening stage, using only the last 5 identified 
frequencies. 
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Fig. 3. Diagrams of the CoRoT light curve of HD 174648, folded 
on the dominant frequencies fi (top) and /2 (bottom), after 
prewhitening all other detected frequencies. The top phase dia- 
gram clearly illustrates the non-sinusoidal shape of the variations 
for /i ; the nonsinusoidal shape for /2 in the bottom diagram is 
less obvious, but is clearly demonstrated by the presence of the 
first harmonic of this frequency at 2/2, labelled /4 in Table[3] 



3.3. Spots and differential rotation 
3.3.1. Observational constraints 

The appearance of two high-amplitude, closely-spaced fre- 
quencies is a natu ral signature of differentially rotating spots 
dLanza et al.ll9 93^, especially if their sec ond harmonic s are also 
prominent, but higher harmonics are not (IClarkell2003h . The fre- 
quency splitting A/ a; 0.05 d ' of HD 174648 suggests that the 
periods of rotation at the spot latitudes differ by ~1%. 

The hypothesis that the photometric variability is due to rota- 
tional modulation implies a high rotational frequency equal to /i 
and /2 at the spots' latitudes. In fact, a confrontation between 
the b reak-up frequen cies for a grid of CL ES dScuflaire et al.l 
20081) models (from iDegroote et al.ll20Toh . shows that /2 ex- 
ceeds the break-up frequency for all models near the location of 
HD 174648 in the Kiel diagram (Fig.|4]i. Nonetheless, if we ac- 
count for the possibility of differential rotation, which is hinted 
at in the spectral profile of Mg n 448 1 A (Sect. l3.1b . then the stel- 
lar surface at higher latitudes could possibly rotate faster without 
exceeding the local break-up velocity. In Fig.m the minimum 
required amount of differential rotation a for a model to reach a 
rotation rate of fi at the pole with Ogq = ^a-it is shown. Here, a 
is defined via 

Q(0) = neq(l - Q-sin-e). (1) 

In the case of polar acceleration a is negative and a minimal 
differential rotation rate is equivalent to a maximal a, thus we 
search for 

_ f^crit ~ fl 
^max — 7T ■ 
^ ^crit 

By placing fi at the pole, we force the differential rotation law 
to be as close to uniform rotation (a = 0) as possible. We see 
that, unless the star is located close to the ZAMS, unrealistically 
high differential rotation rates have to be reached (cf. rates up 
to — 30% reported by Reiners & Schmitt 2003). On the other 
hand, iReiners & Roved (l2004bl) have shown that the differential 
rotation rate increases with effective temperature. The hottest 
stars in their sample have an effective temperature 3 000 K lower 
than HD 174648. Therefore, we restrict the different rotation rate 
to Q- > -40% in the spot modelling in Sect. 13.3.31 

In contrast to the rotation frequency, the projected equatorial 
rotational velocity is dependent on the inclination angle of the 
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star. The only constraint on the spot model imposed by the ob- 
servations is that we cannot have a model with too low a surface 
gravity, to avoid the stellar radius to become too large, and, given 
the observed rotation frequency, the measured Veq sin / to exceed 
the break-up velocity (Fig.|4]i. In the following, we check if all 
data sources can be reconciled with the hypothesis of a differen- 
tially rotating spotted early-type star After defining the general 
model, we will search for and discuss possible solutions. 

3.3.2. The model 

At a given time, the integrated intensity S of the surface of a 
spotted star, with respect to its uns potted intensity ^o, can be 
approximated by jLanza et al ]l200l 

S = SqC Yj A„finh(Mn) [(cs - 1) + Q(Cf + 4//„ - 1)] + So- 

n.fi„>0 

(2) 

The parameter A„ represents the relative surface area of the nth 
spot, which is assumed to be circular, fi„ = cos i/'„ is the co- 
sine of the limb angle of the spot, and C is the unspotted spe- 
cific intensity connected to the limb darkening function h{fi„). 
The parameters Cj and Cf are the starspot and faculae contrast 
parameters, which determine how bright the spot is relative to 
the surroundings. For the faculae, an extra contrast coefficient c'j. 
is introduced because the brightness of faculae is dependent on 
the viewing angle. The ratio between the area of the faculae and 
the starspots is fixed by Q, where Q = means no faculae. We 
chose to fix all spot parameters to the solar values even though 
the physical circumstances in the atmospheres are quite different 
for B stars. The simplicity of our model justifies this choice, as a 
posteriori tests have shown, because the influence of higher con- 
trast parameters, for example, can be mimicked by larger spots. 
In this respect, the physical reality of the formalism we used to 
model the spots is less important than the morphology it intro- 
duces in the light curve. Since our model treats spots as point 
sources, they cannot interfere with each other (i.e. combine or 
deform), and thus larger spots could also be interpreted as darker 
spots. For the Sun, the starspot contrast parameter takes the value 
Cs = 0.67. The facular contrast coefficients are Cf - 1.115 and 

= -0.115, and the ratio between the area of the faculae and 
starspot area is Q - 10. 

In the adopted approach, the limb darkening law is approxi- 
mated by the quadratic law 

/l(jU„) = flo + fliAi + "IfJ^n- 

The integration of Eq.|2]is done over the visible part of the star, 
defined as ju„ > 0, where 

/!„ = COS n sin 0„ + sin ; cos 6„ cos(0„ + f2(0„)f - ^o)- 

The inclination of the star is denoted by ; (by convention, / = is 
pole-on, / = 7r/2 is equator-on), which is the opposite orientation 
of the latitude 0„ (0„ - represents a spot on the equator). Lq is 
the epoch angle where Lq = means that a spot with coordinates 
(9„, (p„) = (0, 0) is at the central meridian of the visible disk when 
t - 0. Finally, the angular velocity Q. is in general dependent on 
the latitude via Eq.[T] 

Spot modelling using only the light curve as a source of in- 
formation is a degenerate problem. It is difficult to apply tradi- 
tional minimization algorithms to the equations to find the cor- 
rect parameters, because there are many local minima in the pa- 
rameter space. Also, it is difficult to find good starting values 



to initialize the fit. In the case of a low number of spots, a ge- 
netic minimization algorithm can be a valuable alternative fit- 
ting method. These kinds of algorithms are generally better at 
exploring large parameter spaces and finding global minima, but 
require more evaluations. To facilitate the search we fix as many 
parameters as possible beforehand (e.g., the spot contrast ratios). 
Nevertheless, the outcome should be regarded as a possible so- 
lution, not a unique one. 

The genetic algorithm to fit Eq.|2] comprises the following 
steps: 

1 . Initialize a sample consisting of thousands of sets of param- 
eters (individuals), each representing a spot model. Choose 
the values of the parameters randomly within appropriate 
search ranges. 

2. Evaluate all spot models, and compute the residuals with re- 
spect to the original light curve. For each solution calculate 
its fitness, i.e. a number representing the quality of the fit. 
If a fit is good, it should have a high fitness. We choose the 
fitness value to be the reciprocal of the;^^- statistic. 

3. Rank the individuals in the sample according to their fitness, 
and select the « best ones. This selection process picks a ran- 
dom fitness threshold between the worst and best individual, 
and retains only those individuals with a fitness higher than 
the threshold. 

4. Keep the best individuals in the sample, and generate off- 
spring to replace the bad individuals. Offspring are gener- 
ated from two good individuals: sometimes switching some 
values, sometimes just copying them (random random-point 
crossover). Also the parents are randomly chosen, meaning 
that some of the good individuals may not generate any off- 
spring, while others have many. Allow a random perturbation 
of all parameters (mutation), but make it small in the begin- 
ning, and slowly raise the perturbation during the evolution. 

Now reiterate steps (l)-(4) to breed successive generations, until 
a stop criterion is reached. This can be either after a certain num- 
ber of generations or at a point where the sample fitness does not 
increase any more. 

3.3.3. Solutions 

The free parameters in the fit describing the star are the equato- 
rial period Peq £ [0.26, 0.43] d and the differential rotation rate 
b - (Qp - Qeq) e [0, 1.5] d"'. Two infinitely lasting spots are 
assumed to be located on the same hemisphere, with longitude 
Ai 6 [0,2;r], latitude 0, e [0,7r/2] and size A,- between 0.001% 
and 0.2% of the surface area. In the genetic fit, the population 
consists of 15 000 individuals, and by default, 100 evolution 
steps are made. On top of the random random-point crossover 
breeding, mutation is added to all breeding phases in the genetic 
fit. This is applied to all parameters aj in the form of normally 
distributed random noise with standard deviation cr - na'^\aj\, 
where ^ - I y. 10"^, a =1.2 and k is the generation number 
This results in a power-law increase of the mutation parameter 
from virtually zero during the first evolution stages, to ~10% at 
the end. The last stages can then be used to check the influence 
of small perturbations on the fit parameters to the fit statistic. 

To investigate the uniqueness and convergence of the ge- 
netic fitting algorithm, we have performed the fit several times. 
From this exercise, we see that many combinations of parame- 
ters give a good fit. For example, with a low differential rotation 
rate, the frequency spacing A/ can be recovered by placing the 
spots at high latitudinal separation. A higher differential rotation 
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Fig. 4. Critical rotation velocities and frequencies, and their relation to stellar parameters for models between 1.5 and 4.1 M©. The 
grey band shows the range of T^ff matching the observed spectra, the + sign shows the location of the best matching model in km s"' 
(see text). Top left panel: critical velocities on the main sequence for CLES models with aov = 0.2. The black line shows the line of 
models matching the observed Vgq sin / with / = 90°. The grey lines show the models reaching their critical rotation frequency for 
different inclination angles for the observed Vgq sin / (/ = 40°, 45°, 60°, 75° and 90°). The colours represent the values of different of 
the following quantities. Top right: /?eq for zero rotation (/?©). Bottom left: Critical rotation frequency (d"' ). Bottom right: Necessary 
differential rotation rate a for a model to reach /i at the pole with Qgq = Odit. 



rate forces the separation to be smaller. This does not mean that 
we cannot formulate any constraints. Indeed, we have to cross- 
correlate the solutions from the fit with the grid of stellar mod- 
els, the measured projected equatorial rotation velocity and the 
spectral profiles. Moreover, the diff'erential rotation rate has to 
be realistic. 

To incorporate all collected restrictions on the star, we de- 
vised the following scheme: 



cannot exceed 20kms^'). We find that all valid models are 
located near the ZAMS. 
3. Finally, we compute synthetic spectra for a selection of ~50 
different valid models, varying in surface gravity, effective 
temperature, inclination angle and rotation frequency. Here, 
we take into account surface distortion using the rigidly 
rotating Roche approximation (ICranmer & Owockil Il995h 



1 . For every model in the grid of evolutionary tracks we know 
the mass and radius in absence of rotation. Thus, we can 
compute the critical equatorial rotation frequency Ociit (uni- 
form rotation approximation) and compare it with the ro- 
tation frequency flgq = l/^'eq retrieved from each separate 
spot model fitted with a realistic differential rotation rate 
(a > -40%). If Qciit < fieq, we consider the fit to be in- 
compatible with the stellar model under consideration. 

2. Next, we can derive the equatorial radius and check if it is 
compatible with the measured Veq sin / and the inclination an- 
gle from the spot model fit (i.e., we require that the difference 



R{9) - ^ 



3R 



p 



Qcosfl 



n + arccos(n cos ff) ] 



(3) 



but we add a differentially rotating velocity field in compli- 
ance with the parameter resulting from the spot model fit. 
We note that this is not entirely consistent with the Roche 
approximation, but the additional surface distortion eff'ects 
due to differential rotation are small. After discretising the 
surface of the star in 400 x 200 surface elements in lon- 
gitude, and latitude, respectively, we compute the synthetic 
spectra in the following way: First, the vectors normal to the 
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surface defined in Eq.|3]are projected onto the line-of-sight. 
Then, a velocity field is constructed or derived and the mag- 
nitude of the velocity field in the line-of-sight is calculated. 
For each surface element, the local gravity log g{0) and tem- 
perature T(6) are derived and the corresponding line profile 
is taken from the library of line profiles. The total line pro- 
file J is then calculated by summing all shifted line profiles 
/ for each surface element, weighted with the amount of flux 
emitted in the Une-of-sight A . The flux is calculated accord- 
ing to the Ivon Zeioell (11924 theorem, and corrected for the 
element area and limb darkening h^, which is derived from 
a library of specific intensities, using again the local gravity 
and temperature, 

J^Yj [^(^' ^^'^'i^^' Si)Iiv, Tu \og gi)] . 

e,(p 

(4) 

To compare the synthetic spectra with the observed spectra, 
we applied the same fitting algorithm as described in Sect. l3.1l 
(Fig.|2] red line). These line profiles are calculated by numer- 
ically summing the contributions of each surface element. For 
each surface element, the spectral features and intensity are de- 
rived from the local effective temperature and local surface grav- 
ity, thus taking gravity darkening effects into account. We note 
that the line profiles can differ more from one another than what 
would be expected from the difference in polar effective temper- 
ature alone. The rotation can lower the effective temperature at 
the equator by several thousand Kelvin, which can have a ma- 
jor effect on the shape of the line, even if the polar tempera- 
ture changes only slightly. The parameters of the best matching 
model are listed in Table|4] and the resulting light curve model 
is compared to the observations in Fig.|5] The rotation law is il- 
lustrated in Fig.|6] and the distribution of the flux, temperature, 
gravity and velocity on the surface are shown in Fig.|7] 

We compare the fit quality of the spotted model with the pul- 
sation model via a reduced statistic, 

2 _ 1 V"' iFmod,A - Fobsj)^ 
^red - Zj ^ • 

Here, cr is determined as V2 times the standard deviation of the 
differenced residual light curve after subtracting the pulsation 
model resulting from a full frequency analysis. Considering only 
two frequencies (7 free parameters), we find xl = 9.14. Adding 
one harmonic of the main frequencies, we end up with;\f2 - 1 -26 
(11 free parameters). The spot model results in xl = 2.61 (10 
free parameters). We conclude that, purely on the basis of fit 
quality, neither the spot hypothesis nor the pulsation hypothesis 
can be excluded. 

Next, we checked if the location of the spots and the star's in- 
clination angle resulting from the best fitting model are compat- 
ible with the amplitude of the variability. Because of CoRoT's 
single-colour white light photometry, we cannot determine the 
spot temperature, and because of the disc integration, we can- 
not infer the morphology of the spots. For simplicity we there- 
fore assumed the spots to be perfectly circular and totally black 
(i.e., blocking all radiation and not experiencing any faculae). 
We then placed these dark patches on the distorted surface of 
the model star, and computed the minimum size of the spots to 
induce ~0.18% and ~0.13% flux variations. The largest spot at 
86.6° is compatible with a spot radius of ~4° (parameter ^'j in 
Table|4]i, or a surface coverage of ~Q.1%. The smallest spot at 
75.9° should have a radius of ~3° corresponding to a surface 



Table 4. Spot model parameters for HD 174648. 



Parameter 


Value 


M(Mo) 


2.4 ±0.1 




1.7 


reff.polc (K) 


10 950 + 500 


loggpole (dex) 


4.37 ± 0.2 


P(d) 


0.34 ± 0.02 


a 


-0.32 ± 0.03 




0.87 


^pole (^crit) 


1.14 


Veq sin I (km s"') 


283 


in 


76.3 ±4 


01 n 


75.9 ±4 


62(1 


86.6 ± 4 


Si (surface fraction) 


0.0135 ±0.005 


$2 (surface fraction) 


0.0385 ± 0.003 


s[ (surface fraction) 


0.0010 {r ^ 4°) 


s'y (surface fraction) 


0.0049 (r ^ 3°) 




2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 
Rotation frequency (d""^ ) 



Fig. 6. Differential rotation law for the best fitting spot model 
(solid line). In dashed lines, the breakup rotation frequency is 
shown as a function of latitude. The locations of the two spots 
are shown as white circles. 



coverage of ~0.05% (parameter s'2 in Table|4|l. These are plausi- 
ble values; the spots can be small compared to the surface area, 
and do not need to overlap. In fact, a reduction of 10% of the 
original flux at the spot locations can already induce the same 
variability with only a marginal increase of the spot size. These 
values are smaller than those resulting from the spot model fit 
(TableHl parameters S[ and ^2), where we assumed solar values 
for all spot parameters. 

3.3.4. Discussion 

The two dominant frequencies in the CoRoT light curve of 
HD 174648 both have a clear second harmonic, which is a typ- 
ical signature of spot-like features. The fact that we do not see 
a constant part in the phase shapes of the light curve folded on 
the dominant frequencies implies that the spots never fully dis- 
appear behind the stellar disk. This means that we either see the 
star nearly equator-on with two circumpolar spots, or that the 
star has a significant inclination with respect to our line of sight 
and the spots can be anywhere on the visible hemisphere. The 
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3190 
Time (CJD) 



Fig. 5. Spot model fit and residuals. Top panel: Black dots and lines show the observed CoRoT light curve, corrected for disconti- 
nuities. In grey, the best fitting spot model with parameters listed in Table|4]is shown. The inset is a zoom on the middle part of the 
light curve. Bottom panel: Residuals after subtracting the spot model. 




Fig. 7. Stellar parameters across the rotationally distorted sur- 
face of the star, and location of the spots. The paths of the spots 
are shown in black lines. The equator of the star is shown as a 
white line. Top left: emitted flux (linear scale). Top right: rotation 
velocity (km s"'). Bottom left: effective temperature (K). Bottom 
right: logarithmic surface gravity (cgs) (coloured version only 
online). 



first possibility is favoured by the high projected rotational ve- 
locity derived from the Mgii4481A line profile. The hypothe- 
sis of circumpolar spots is also supported by the rotation fre- 
quency at the spot latitudes. The frequencies fx and /2 are too 
high to be equatorial rotation frequencies, since they are higher 
than the break-up frequency for any realistic stellar model com- 
patible with the fundamental parameters of the star. Fast rotating 



circumpolar spots are then compatible with an equatorial rota- 
tion frequency below break-up velocity and a high differential 
rotation rate with polar acceleration. A sign of high differen- 
tial rotation with polar acceleration is also found in the shape 
of the M g II 4481 A line profile. According to lReiners & Schrniti 
( l2002al) . polar spots can mimic differential rotation with polar 
acceleration in spectroscopic line profiles, but this degeneracy 
is lifted based on the light curve morphology. Uniform rotation 
would remove the beating pattern since the two spots would have 
the same rotation period. 

We have conceptually proven the possibility of the two-spot 
model to explain the light curve of HD 174648, and we have 
shown that the resulting parameters are compatible with mod- 
els of stars close to the ZAMS. The solution is in agreement 
with the observed spectral lines. Four key features have to be 
explained to make the spot hypothesis plausible from a purely 
physical vantage point: rapid rotation, a high differential rotation 
rate, polar acceleration and the existence of circumpolar spots. 
Rapid rotation is most prominent in stars of spectral type A and 
B, of which the former have the highest fractional rotation ve- 
locities on average, and the latter the highest absolute rotation 
velocities. The polar acceleration is compatible with r apid rota- 
tion in stars without a significant convective outer layer dClemenj 
[19691) : convection with viscous braking implies equatorial accel- 
eration (e.g., the Sun), but in the absence of significant convec- 
tive layers, radiative braking takes over as the most important 
braking mechanism, and results in polar acceleration. Due to the 
high effective temperature of HD 174648, we expect this star 
only to have a shallow convective layer, if any, an d thus rapid 
rotatio n and polar acceleration are indeed possible. IZorec et all 
(1201 ih show that rapid rotation can lead to differential rotation in 
early type stars, and Hussain (2 002h andll gik et al. (2007 ) show 
that rapid rotation can lead to a concentration of spots at polar 
latitudes, because of an efficient transport of flux from lower lat- 
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itudes via meridional circulation. The same mechanism can also 
transport the magnetic field to higher latitudes and induce strong 
polar magnetic fields, benefiting large, long-lived spots. Finally, 
iHussainI (l2002h shows that the difl'erential rotation rate increases 
with stellar mass, though we have to be careful not to overin- 
terpret her simulations, since they were carri ed out for cooler 
stars (m asses up to 1.7 M©). On the other hand. iReiners & RoveJ 
(l2004bl) have analysed a sample of A stars, and found evidence 
for a high differential rotation rate of ~3Q% in the earliest-type 
object in their sa mple (HD 60555, A6V, vsin; = 114km s '; 
iRover et al.ll2007l) . 

Nonetheless, some problems remain. If only a very small 
convective region exists, how can a m agnetic field be gene r- 
ated strong enough to create spots (cf. ICantiello et al.ll20Tol) ? 
Even though the above literature studies are supportive of 
our interpretation, a strong correlation between the presence 
of mag netic fields and low rotation rate was previously ob- 
served jBriquet et alj 120071) . which is not in agreement with 
the fast rotation of HD 174648. Magnetic fields and/or spots 
have only been detected in the slowly rotating chemically pecu- 
liar Ap/Bp stars, HgMn st ars (Wolff 1983; Briquet et al. 2001, 
12007'. Alecian et alJ l2009t fBriquet etalJ 1201 01) '. and in SPBs 
(iHubrig et al. 1 120091) . Thus, if a magnetic field is present, it 
should not only be weak, but also fast-evolving and probably 
non-dipolar, to allow for the existence of close circumpolar spots 
and differential (fast) rotation. This scenario is not reconcil- 
able with what has been found for th e magnetic Ap/Bp stars 
dShorlin et al.l2002l:lAuriere et alJ2007h . Yet, rotational modula- 
tion has been proposed for line profile variations in other rapidly 
rotating stars, sometime s in combination with pu lsations (e.g., 
the ySCep binary kSco, lUvtterhoeven et al"1l2005h . Another ex- 
planation for the spots is to invoke inhomogeneities in surface 
chemical abundances. In this case, the problems linked to Bp- 
like magnetic fields vanish, since these inhomogeneities are ob- 
served in stars for which no large-scale or strong magnetic fields 
have been detected (e.g., iWade et al.i,2006.; . Auriere et al. 2010,). 

A final remark concerns the possibility of contamination of 
the CoRoT light curve by variability of a cooler star in the 
same field of view. A foreground spotted K5V star can induce 
an amplitude of 1800ppm in the light curve. At the distance of 
HD 174648 or further away, the spot amplitude has to be larger 
than 0.2 mag to cause the observed variabiUty, which is unlikely. 
The detected frequency, however, would imply an equatorial ro- 
tation velocity of 170 km s"' , which can only be obtained in rare 
cases of a tidally locked compact binary (e.g., a K5 V star around 
a white dwarf with a semi-major axis of ~2Rq). 

4. The pulsating stars 

In terms of light curve morphology, features in the frequency 
spectrum (e.g., the presence of harmonics) and rotation rate 
(veq sin; = 271 ± 16km s-i), HD 170935 resembles HD 174648 
(Fig.|9]i- The surface gravity suggest that it is slightly more 
evolved, frequency values themselves are considerably lower 
(~0.70d ' compared to ~3.8d '), and the number of frequen- 
cies higher with respect to HD 174648. Explaining this light 
curve in terms of spots would require a complex pattern of nu- 
merous spots. Nonlinear gravity mode pulsations are therefore 
a more likely cause of variability, making H D 170935 a fast ro- 
tating and cool SPB star, just as HD 121 190 jAerts & Kolenberd 
120051) . 

The four other cool B stars in the CoRoT sample (Figs lA.U 
IA.7l i contain fewer frequencies than the two stars discussed pre- 
viously. Moreover, at low amplitudes, residual instrumental ef- 
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Fig. 8. Short time Fourier transform of the CoRoT light curve of 
HD 181440 (window length of 30 d). The two highest frequen- 
cies are constant throughout the time series. The low frequency 
could result from rotational modulation, but a beating pattern 
cannot be excluded (amplitudes in ppm). 

fects can be important, diminishing the number of detected fre- 
quencies even further. Except for the frequencies with very high 
S/N, we assume that the bulk of the low frequencies are re- 
lated to instrumental noise. The only star in the remainder of 
the sample with a clear multiperiodic frequency spectrum is 
HD 181440 (B9III). The high S/N frequencies and candidate 
combination frequencies are listed in Table|5] The presence of 
a multiplet around /i = 0.44414 + 0.00014d"' (containing /4,/5 
and /e) hints for existence of variable amplitudes, frequencies, 
or phases, but a beating pattern cannot be ruled out. No signif- 
icant harmonics are extracted. The possibility of a multiplet is 
excluded from the low value of the spacing and the projected 
rotational velocity. The frequencies /2 and fj, are isolated and 
therefore assumed to be pulsations. 

In the frequency spectra of the two slowly rotating stars 
HD 182198 (B9V) and HD 49677 (B9), only one frequency is 
clearly present. The frequencies / - 0.13976 + 0.00001 d"' 
(in HD 182198) and / = 1.3364 + 0.0009d-i (HD49677), are 
good candidates for an SPB gravity mode and could therefore 
be probes for the borders of the SPB instability strip. The oc- 
currence of monoperiodic SPBs, or of such pulsator s with a 
highly dominant mode, is common (iDe Cat et al.ll2005h . Finally, 
no clear stellar frequencies were recovered from the moderate 
rotator HD 46179. 

5. Conclusions 

The cool side of the theoretical SPB instability strip harbours 
many different types of stars, some of them exhibiting clear vari- 
ability. SPB stars co-exist in this part of the HR diagram with 
Be stars, Bn stars. Bp stars and HgMn stars. The reasons for 
the differences among these classes of stars are not well under- 
stood, but rotation and magnetic fields must play important roles. 
We have argued that pulsations need not be the (only) source 
of variability for all these stars. Instead, rotational modulation 
is an alternative, at least for the fast rotators. For slow rotators 
or objects with low surface gravities, the rotation periods are of 
the same order as the expected pulsation periods, thus confusion 
occurs if only white-light photometric time series are available. 
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Table 5. Frequency analysis of the CoRoT light curve of HD 181440, after correcting for discontinuities. Only the frequencies with 
a 5 /A^ > 20 are shown, and candidate combination frequencies. 



ID 


f 


1 


a 
ppm 


<P 

rad/27r 


S/N 




/i 


0.444136 ± 


0.000014 


72.3 ± 0.3 


0.432 ± 0.004 


39.58 




/2 


0.620652 ± 


0.000013 


69.6 ± 0.3 


-0.332 ± 0.004 


40.50 




h 


1.166012 ± 


0.000017 


47.0 ± 0.2 


0.441 ± 0.005 


30.60 




/4 


0.449199 ± 


0.000024 


33.7 ±0.2 


-0.354 ± 0.007 


22.91 




h 


0.430491 ± 


0.000026 


27.9 ± 0.2 


0.445 ± 0.008 


21.45 




h 


0.410365 ± 


0.000026 


26.9 ± 0.2 


0.053 ± 0.007 


22.11 




fi 


0.543475 ± 


0.000034 


20.2 ± 0.2 


-0.241 ± 0.010 


16.89 


/1-/2 + 0.001890 


h 


0.180443 ± 


0.000055 


12.2 ±0.2 


-0.189 ±0.016 


11.66 


/2 - /i + 0.003927 


f9 


0.174743 ± 


0.000097 


7.6 ±0.2 


0.219 ±0.028 


7.20 


/2 -/4 + 0.003291 
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Fig. 9. Excerpt of the CoRoT light curve of HD 170935, showing similar variability as HD 174648 (grey line shows a fit to the light 
curve, black dots represent to CoRoT data). 
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Fig. 10. Frequency spectrum of the CoRoT light curve of HD 170935 above 0.1 d ' (black). The significant frequencies and fitted 
amplitudes above 10 ppm are indicated in grey. 



Although time-frequency analyses hold the potential to distin- 
guish between the two phenomena, the frequency spectrum of 
gravity modes is expected to be dense, thus long term continu- 
ous observations are needed to disentangle beating patterns from 
time-dependent behaviour. 



Despite the precision and duty cycle of single-band space- 
based photometry, we have shown that a unique and consistent 
description of late B star variability is diflicult to obtain solely 
based on this source of data, even when multicolour photome- 
try or spectroscopy are available to constrain the fundamental 
parameters such as the effective temperature and surface gravity 
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of the star. We have also seen that a good determination of the 
rotation frequency is important. Additionally, multicolour pho- 
tometric time series or spectroscopic time series are necessary to 
determine the origin of variability of these cool B stars, as spots 
and pulsations leave distinct signatures in the amplitude ratios 
for different photometric bandpasses, and in line profiles. 

We found evidence for the existence of differential rotation at 
the surface of one fast rotator of late spectral type B. This holds 
information on the internal rotation profile if we can pinpoint 
the cause of the nonrigidity, and can thus be a tool to probe the 
interior of massive stars. Only cool B stars were studied, but 
there is no reason to assume similar stars cannot exist with hotter 
temperatures. In fact, we could hope for the existence of SPB 
stars with inhomogeneous surface structures, which would give 
independent constraints on the rotational behaviour of these stars 
to compare with those obtained from asteroseismology. 

The origin of the spots remains an unanswered question. A 
strong. Bp-like magnetic field is incompatible with fast rotation 
and the need for a non-dipolar field to e xplain the close cir - 
cumpolar spots to be of magnetic origin (lAuriere et al ] 120071) . 
On the other hand, a weak magnetic field and the existence of 
differential rotation could more easily explain the fast rotation, 
but would invoke the need for fast evolving magnetic fields to 
generate the spots. Chemical spots are a valuable alternative, 
since they appear not to need the presence of magnetic fields, 
though there is no consensus in the literat ure on an alterna- 
tive underlying mechanism to generate them dHubrig & CastelUl 
l200ltlWade etaHIIOOeHAuriere et al.ll2010h . 
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Appendix A: Light curves and frequency spectra of 
B8/B9 stars observed witli CoRoT 
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Fig. A.l. Excerpt of the CoRoT light curve of HD 1 8 1440, showing multiperiodic variability (grey Une shows a fit to the light curve, 
black dots represent to CoRoT data). 
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Fig. A.2. Frequency spectrum of the CoRoT light curve of HD 181440 above 0.1 d ' (black). The significant frequencies and fitted 
amplitudes above 10 ppm are indicated in grey. 
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Fig. A.3. Excerpt of the CoRoT light curve of HD 182198, with only one clear clear frequency (grey line shows a fit to the light 
curve, black dots represent to CoRoT data). 
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Fig. A.4. Frequency spectrum of the CoRoT light curve of HD 182198 above 0.1 d ' (black). The significant frequencies and fitted 
amplitudes above 10 ppm are indicated in grey. 
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Fig. A.5. Excerpt of the CoRoT light curve of HD 49677, with only one clear frequency detected (grey line shows a fit to the light 
curve, black dots represent to CoRoT data). 
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Fig. A.6. Frequency spectrum of the CoRoT light curve of HD 49677 above 0. 1 d ' (black). The significant frequencies with ampli- 
tudes above 10 ppm are indicated in grey. 
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Fig. A.7. Excerpt of the CoRoT light curve of HD 46179. No clear signs of variability are detected (grey line shows a fit to the light 
curve, black dots represent to CoRoT data). 
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Fig. A.8. Frequency spectrum of the CoRoT light curve of HD 46179 above 0.1 d ' (black). The significant frequencies with ampli- 
tudes above 10 ppm are indicated in grey. 



